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Abstract. We analyzed the spectroscopic and photometric database of the 5th data release of the Sloan Digital Sky Survey 
(SDSS) to search for evidence of the quasar reddening produced by dust embedded in intervening damped Ly a (DLA) systems. 
From a list of 5 164 quasars in the interval of emission redshift 2.25 < z e < 3.5 and SDSS spectra with signal-to-noise ratio SNR 
> 4, we built up an "absorption sample" of 248 QSOs with a single DLA system in the interval of absorption redshift 2.2 < 
z a ^ 3.5 and a "pool" of 1959 control QSOs without DLA systems or strong metal systems. For each QSO of the absorption 
sample we extracted from the pool a subset of control QSOs that are closest in redshift and magnitude. The mean color of 
■ tms subset was used as a zero point to measure the "deviation from the mean color" of individual DLA-QSOs, A,. The colors 

' were measured using "BEST" ugriz SDSS imaging data. The mean color excess of the absorption sample, (£), was estimated 

by averaging the individual color deviations A,. We find (E(r - z)) = 27 (±9) x 10~ 3 mag and {E{g - z)) = 54 (±12) x 1CL 3 
mag. These values are representative of the reddening of DLA systems at z a ~ 2.7 in SDSS QSOs with limiting magnitude 
r - 20.2. The detection of the mean reddening is confirmed by several statistical tests. Analysis of the results suggests an origin 
of the reddening in dust embedded in the DLA systems, with an SMC-type extinction curve. By converting the reddening into 
rest-frame extinction, we derive a mean dust-to-gas ratio (A v /N(Hl)) « 2 to 4 x 10~ 23 mag cm 2 . This value is ^ -1.25 dex 
7—i ' lower than the mean dust-to-gas ratio of the Milky Way, in line with the lower level of metallicity in the present DLA sample. 

: 
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1 . Introduction At redshift z < 1.6 about 20 DLA host galaxies have been 

imaged in the quasar field, showing a variety of morphological 

Quasar absorption-line systems allow us to probe the physi- types (Le Brun et al. 1997, Rao et al. 2003, Chen et al. 2004). 

cal properties of intergalactic/interstellar gas over a large frac- At higher re dshift, where the bulk of known DLA systems are 

tion of the Hubble time, back to the epoch of quasar forma- i oca ted, the imaging identification of the galaxy is quite dif- 

tion. The quasar absorbers with neutral hydrogen column den- ficult and most of me studies on the nature of DL A galaxies 

sity N(H i) > 10 2() 3 atoms cm 2 are believed to originate in m based on spectroscopic observations. In spite of a long se- 

the interstellar gas of intervening galaxies (Wolfe, Gawiser, ries of studies on their kin ematics (Prochaska & Wolfe 1997, 

& Prochaska 2005), rather than in the intergalactic medium. 1998)i chemical abundances (Pettini et al. 1994; Lu et al. 1996; 

The presence of a strong Lya profile with "damping" wings p roc h a ska & Wolfe 1999; Molaro et al. 2000; Dessauges- 

extending well beyond the Doppler core of the line is a dis- Zavadsky et al. 2006), and physical properties (Petitjean et 

tinctive feature that gives the name to this class of absorbers al 2 000; Wolfe, Prochaska & Gawiser 2003), the exact na- 

(Wolfe et al. 1986). Spectroscopic studies of "damped Lya" ture D f DLA galaxies is still open to debate (Wolfe & Chen 

(DLA) systems yield unique and accurate information on gas- 2 006). Studies of DLA systems allow us to probe faint galax- 

rich, high-redshift galaxies observed in absorption (hereafter ies that happen to lie in the direction of the quasar. It is not 

DLA galaxies ). clear whether there is a continuity or not between the proper- 
ties of absorption-selected DLA galaxies and photometrically 
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selected Ly-break galaxies, which are biased in luminosity (e.g. 
M0ller et al. 2004). 

In the present work we focus our attention on the dust 
component of DLA systems. By analogy with studies of local 
galaxies, we expect interstellar dust to be a pervasive compo- 
nent of DLA galaxies. The properties of dust are poorly un- 
derstood, but may play a critical role in several investigations 
aimed at understanding the nature of DLA galaxies. In fact, 
dust is expected to play a major role in a variety of physical 
processes and to affect the measurement of observational quan- 
tities. 

An example of the importance of dust in governing the 
physical processes at work in DLA galaxies comes from the 
study of the C n* 1335.7 A line, which can be used to estimate 
the star formation rate assuming that the heating is dominated 
by photoelectric emission from dust grains (Wolfe, Prochaska 
& Gawiser 2003). The efficiency of the heating mechanism and 
the accuracy of the SFR determined in this way depend on the 
dust-to-gas ratio, which is poorly constrained. Also the pro- 
duction of molecular gas is believed to be influenced by the 
presence of dust, which acts as a catalyst of H2 formation. 

An example of the impact of dust on the observations of 
DLA systems comes from the different estimates of chemi- 
cal abundances obtained from refractory and volatile elements 
(Pettini et al. 1994; Vladilo 1998, 2002; Centurion et al. 2000; 
Hou et al. 2001, Prochaska & Wolfe 2002). A significant frac- 
tion of refractory elements is depleted from the gas phase, 
where they can be detected via absorption line spectroscopy, 
to the dust component, where they escape detection with this 
technique. Studies of galactic chemical evolution are based on 
measurements of abundance ratios which may be affected by 
differential dust depletion. The inferences of galactic chemical 
evolution models on the nature of DLA galaxies depend on the 
exact amount of dust depletion in these systems (e.g. Calura et 
al. 2003, Dessauges-Zavadsky et al. 2006). 

Another observational effect that we expect to occur is the 
absorption and scattering of the photons of the quasar contin- 
uum by dust embedded in the intervening DLA system. This 
extinction process is more efficient at shorter wavelengths and 
is expected to redden the quasar colors. In the most extreme 
cases the extinction could lead to the obscuration of the back- 
ground quasars (Ostriker & Heisler 1984; Fall & Pei 1989, 
1993). In the more general case, the extinction may induce a 
selection bias acting against the detection of metal-rich galac- 
tic regions in magnitude-limited surveys of quasars (Prantzos 
& Boissier 2000, Vladilo & Peroux 2005). Studies of radio- 
selected quasars surveys suggest that the impact of this effect 
on the statistical properties of DLA systems is modest (Ellison 
et al. 2001, Akerman et al. 2005, Jorgenson et al. 2006), but the 
size of these surveys is not sufficiently large to reach firm con- 
clusions. The existence of empty fields without optical identi- 
fications in the radio-selected survey of Jorgenson et al. (2006) 
suggests that a fraction of quasars may indeed be obscured. A 
better understanding of the amount and type of dust present in 
DLA systems is fundamental to assess the importance of the 
dust extinction effect. 

Because of the potential impact of dust in studies of DLA 
systems, a significant observational effort has been dedicated 



to prove its existence and to understand its properties. At red- 
shift z < 2 definitive evidence for DLA dust has been found in 

one line of sight where the dust extinction bump at 2175 A and 
the silicate absorption at 9.7 fi have been detected (Junkkarinen 
et al. 2004, Kulkarni et al. 2007). Evidence for reddening due 
to DLA absorbers at redshift z < 2 has been found in a few 
quasars with metal rich DLA systems (Vladilo et al. 2006). At 
z > 2 the evidence for dust in DLA systems is mostly based 
on studies of elemental depletions (Pettini et al. 1997, 2000; 
Hou et al. 2001; Prochaska & Wolfe 2002; Vladilo 1998, 2004; 
Dessauges-Zavadsky et al. 2006) and on the existence of gen- 
eral trends between depletion, metallicity and H2 molecular 
fraction (Ledoux et al. 2003, Petitjean et al. 2006). In order to 
establish more firmly the presence of dust at z > 2 and to under- 
stand its properties, it is fundamental to complement the studies 
of depletion with measurements of quasar reddening. A claim 
of reddening detection was reported by Pei et al. (1991), based 
on the analysis of a sample of 13 DLA-quasars. This claim was 
not confirmed by subsequent work, including a photometric 
study of the colors of radio-selected quasars with and with- 
out DLA systems (Ellison et al. 2005) and a study of quasar 
spectra of the 2nd data release of the Sloan Digital Sky Survey 
(SDSS) (Murphy & Liske 2004). The remarkably low upper 
limit of reddening found by Murphy & Liske, (E(B-V)) < 0.01 
mag in the absorber frame, indicates how challeging this type 
of measurement is. On the other hand, the detection of quasar 
reddening due to Can systems (Wild & Hewett 2005, Wild et 
al. 2006) and Mgn systems (York et al. 2006, Menard et al. 
2007) suggests that also DLA systems should produce a red- 
dening signal. 

The aim of the present work is to exploit the large database 
of the 5th SDSS Data Release, in conjunction with novel fea- 
tures in the analysis, to detect the mean reddening of DLA- 
QSOs at z > 2. As in previous work we use the spectroscopic 
database in the process of selection of the quasars of the ab- 
sorption and control samples. At variance with previous work 
on DLA reddening, the quasar colors are measured making use 
of the photometric database. The photometric measurement of 
the reddening avoids some critical steps inherent to the spectro- 
scopic method, such as the photometric calibration of the spec- 
tra and the tracing of the quasar continuum. Uncertainties in 
the photometric calibration, even if small (Adelman-McCarthy 
et al. 2007), translate into uncertainties in the slope of the con- 
tinuum. 

In Section 2 we describe the process of selection of the 
quasar samples. Novel features of our approach include (i) the 
use of spectra with the same signal-to-noise ratio for the selec- 
tion of DLA quasars and control quasars, (ii) the rejection of 
low-redshift absorption systems that may affect the measure- 
ment of the reddening at high-redshift, and (iii) the rejection of 
quasars with multiple DLA systems. Thanks to these last crite- 
ria we are in the position to search for correlations between the 
color excess and the properties of individual DLA systems. In 
Section 3 we explain the method adopted to measure the mean 
reddening of the absorption sample. The measurement is pre- 
sented in Section 4 and its interpretation is discussed in Section 
5. 
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2. The quasar samples 

The starting list of quasars was extracted with the SDSS DR5 
spectroscopic query format. A total of 7294 objects with spec- 
troscopic constraints "QSO" or "HiZ_QSO", emission redshift 
2.25 < z e < 3.5 and imaging constraints "Point Sources" were 
extracted from the query. Quasars at lower redshift do not yield 
sufficient redshift path to search for DLA absorptions. The up- 
per limit on z e comes from the requirement that the photometric 
band r falls redwards of the quasar Ly a emission (see Section 
3.2). 

The detectability of the spectral signatures suggestive of 
reddening depends on the signal-to-noise ratio (SNR) of their 
SDSS spectra. In order to perform a homogeneous comparison 
between absorption and control samples we defined a common 
threshold value of SNR. In practice, we adopted the criterion 
SNR > 4, where SNR is the mean signal-to-noise ratio per pixel 
calculated in the spectral window 1440-1490A in the rest frame 
of the quasar. We found 5164 quasars with 2.25 < z e < 3.5 and 
SNR > 4 in the adopted spectral window. 

From this list of quasars we extracted the "absorption sam- 
ple" and the "control sample". The absorption sample is the 
subset of quasars with a single DLA system and without any 
other potential sources of foreground reddening. The control 
sample is the subset of quasars without DLA systems or any 
other potential sources of foreground reddening. 

The identification of DLA systems and potential sources of 
reddening is based on the analysis of the quasar spectra, as ex- 
plained below. Concerning the DLA identification, we adopt 
tight criteria in order to select only genuine DLA systems. 
Concerning the other sources of reddening, we instead adopt 
broader criteria since we want to be sure to reject any potential 
source that may affect the reddening measurement. 

The list of signatures of potential sources of reddening that 
we considered includes: strong HI absorptions at redshift z + 
Zdla, strong, low-ionization metal systems at z + Zdla, Broad 
Absorption Line (BAL) features, and strong self-absorptions in 
the quasar Ly a emission. Strong metal lines of low ionization 
allow us to search for low-redshift intervening systems not di- 
rectly detectable as a Ly or line because of the limited spectral 
coverage of SDSS spectra. For instance, Mgll absorbers are de- 
tectable down to redshift z = 0.95 in quasars at z e = 3.5, and 
down to z — 0.4 in quasar at z e = 2.25. Other low-ionization 
metal lines redwards of the Ly a emission (e.g. Sin 1526 A, 
Fe n 2600 A) can trace potential sources of reddening at low 
redshift. BAL features and strong Ly a self-absorptions arise 
in the quasar environment and represent a reason of concern in 
the present analysis since they may be associated with dust. In 
fact, BAL quasars tend to show peculiar colors relative to non- 
BAL ones (see e.g. Trump et al. 2006), possibly as the result of 
dust present in the broad line region. Strong HI self-absorptions 
at z — z e are suggestive of large amount of neutral gas, poten- 
tially associated with dust. 

2.1. The selection process 

To build up the absorption sample we started by identifying 
Damped Ly a systems with the automated algorithm described 



by Prochaska et al. (2005; hereafter PHW05) applied here to 
the DR5. Only systems with Ly a absorption blueshifted by at 
least 3000 km s _1 relative to the quasar Lya emission^ were 
considered. Quasars with BAL flags=l and 2 according to the 
criteria explained in PHW05 were rejected. In this way we ob- 
tained a lisi of 422 DLA systems in quasars with z e < 3.5 and 
SNR > 4 in the rest-frame window 1440-1490A. 

Starting from this list the final absorption sample was built 
up in the following way. We first discarded 69 quasars with 
strong Mgll absorbers at redshift zwisii + Zdla- In practice, due 
to the limitation of the SDSS wavelength coverage, we were 
able to search for absorbers with ZMgii < 2.2. To implement this 
rejection process we used the SDSS catalog of strong Mgll ab- 
sorbers with rest-frame equivalent width > 1 A (Prochter et al. 
2006) updated to DR5 (Prochter 2007, priv. comm.). We then 
discarded 54 quasars with multiple DLA systems. Finally we 
rejected 30 quasars with a strong HI absorption at z + Zdla- 
For this purpouse we used a list of candidate HI absorptions 
automatically generated. We rejected cases with /V(HI) > 10 20 
atoms cirT 2 within the errors. 

After this automatic process we visually inspected the re- 
maining 269 quasar spectra to search for spectroscopic fea- 
tures undetected by the automated algorithms. A dozen of 
quasars with strong Mgn absorption not associated with the 
DLA systems were discarded in this way, as well as a few cases 
with anomalous spectra (e.g. extremely narrow Ly a emission). 
Other unidentified weak features were tentatively attributed to 
metal lines associated with the DLA system. 

At the end of the rejection process we obtained a list of 248 
DLA-quasar pairs representing our absorption sample, listed in 
Table 1. At the mean emission redshift of this sample, z e — 3, 
the quasars are free of foreground (Mgll) systems down to ab- 
sorption redshift z a ~ 0.75. The chance of intersecting an addi- 
tional DLA system at lower redshifts is low, so that the result- 
ing DLA-quasar list is a good approximation of an ideal sample 
of quasars with a single intervening DLA system at high red- 
shift (2.2 < z a < 3.5). 

To construct the control sample we started again from the 
list of 5164 DR5 quasars with 2.25 < z e < 3.5 and SNR > 4 
in the spectral window 1440-1490A. We then rejected quasars 
with spectral signatures of potential sources of reddening, us- 
ing the same automated algorithms adopted for the absorp- 
tion sample. In this way we discarded 270 BAL quasars (BAL 
flag=l or 2), 997 quasars with strong Mgn lines (rest frame 
equivalent width > 1 A), and 669 quasars with strong H i sys- 
tems (including DLA systems closer than 3000 km s _1 to the 
Ly a emission). 

The remaining 3228 quasars show, in many cases, absorp- 
tion systems undetected by the automated algorithms. This is 
due to the fact that the algorithms do not search for absorption 



1 This criterion minimizes the likelihood that the DLA absorption 
is physically associated with the quasar. A study of the DLA systems 
closer than 3000 km s _1 to the Ly a emission is presented in a separate 
work (Prochaska, Hennawi & Herbert-Fort 2007). 

2 The complete list of DLA systems of the DR5 includ- 
ing quasars with z e > 3.5 an be found at the web page 
http://www.ucolick.org/~xavier/SDSSDLA 
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Table 1. Data for the quasars of the absorption sample and measurements of the deviations from the mean color A(r— z) and A(g- 
z). Measurements in the (g - z) color index are only given for quasars with z e < 3.4 (see Section 3.2.2). The first 20 DLA-quasars 
are shown. The full table is available in electronic form (see also web page http://adlibitum.oats.inaf.it/vladilo/tables.html). 
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" Width of the bin in emission redshift adopted for selecting control quasars. 
* Width of the bin in z magnitude adopted for selecting control quasars. 

' Mean color and dispersion of the control quasars selected at the same redshift and magnitude of the DLA-quasar. 
d Deviation from the mean color of the DLA-quasar. 



lines in spectral regions with insufficient SNR (e.g. portions 
of the Lyman a forest particularly crowded) or close to strong 
quasar emissions (e.g. the O vi emission). While this limitation 
is not critical for the selection of the absorption sample (only 
bona fide Damped systems are automatically identified), it is a 
reason of concern for the selection of control quasars free of 
any potential source of reddening. We therefore proceeded to 
a visual classification of the quasar spectra in a search for in- 
tervening absorbers not detected automatically. Metal lines and 
strong H i systems were tentatively identified by visual inspec- 
tion with the following criteria. 

Narrow absorption lines redward of the QSO Lya emission 
with central depth > 0.2 relative to the adjacent continuum and 
located away from the CIV/SilV quasar profiles were tenta- 
tively attributed to low-ionization metal absorptions. The SNR 
preselection was in general sufficient to identify features with 
central depth > 0.2 in the full spectral range. For the typical 
FWHM that can be appreciated by eye in the SDSS spectra 
( > 20 A) a central depth of 0.2 corresponds to an equivalent 

width of 4 A in the observer rest frame. This corresponds to 
a rest-frame equivalent width of a; 1 A (« 2 A) for absorbers at 
redshift z a = 3 (z a = 1). Metal lines with equivalent widths of 
this magnitude are generally saturated and suggestive of strong 
metal absorptions potentially associated with DLA systems. 



Saturated absorptions in the Ly a forest with FWHM > 40 

A and zero residual intensity were attributed to strong HI sys- 
tems. For the typical absorption redshift detectable in our spec- 
tra, this threshold roughly corresponds to an equivalent width 
of « 10 A in the rest frame of the QSO. This value is tradition- 
ally used to identify candidate Damped Lya absorptions in low 
resolution spectra (Wolfe et al. 1986). For features lying on the 
wings of strong quasar emissions the criterion was relaxed to 
FWHM > 20 A, given the fact that part of the absorption is 
probably washed out by the emission itself. 

In the course of the visual inspection process we also 
searched for quasars with strong Ly a self-absorption and BAL 
quasars not identified automatically. 

As a result of the visual inspection process we rejected 37 
quasars with BAL features, 1074 suspect metal system^, 141 
potential DLA systems (including absorptions in the proxim- 
ity of the Ly a emission), and 17 cases with a strong Ly a self- 
absorption. In this way we were left with a sample of 1959 con- 
trol quasars listed in Table 2. This is a reasonable approxima- 
tion of an ideal sample of quasars without intervening sources 
of reddening. As we explain in the next section, this sample 
is used as a "pool" from which we extract subsets of control 

3 This number is consistent with estimates of the number of inter- 
vening DLA systems at z„ == 2.2 not detectable in Ly a, but potentially 
detectable as metal systems. 
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quasars specifically designed for each quasar of the absorption 
sample. 

Table 2. Pool of control quasars used in the present in- 
vestigation (see Section 2). The first 20 quasars are shown. 
The full table is available in electronic form (see also 



http://adlibitum.oats.inaf.it/vladilo/tables.htmI). 
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Fig. 1. color excess E(r - z) in the observer frame calculated 
for an absorber at redshift z a with rest-frame extinction A y = 1 
mag. The lines correspond to an SMC-type extinction curve 
(magenta, solid line; Pei 1992, Gordon et al. 2003) and three 
MW-type extinction curves (red lines; Cardelli et al. 1998) 
computed at three values of R = Ay/E(B - V) (dots: R = 2.5; 
short dashes: R = 3.1; long dashes: R = 4.0). 



" Mean signal-to-noise ratio in the rest-frame spectral window 1440- 
1490A. 



where (A(y - x)) is the average of the «,■ values A,-. The validity 
of Eq. (O is discussed in Appendix A. 



3. The method 

3.1. Derivation of the mean color excess 

Given two measurements m x and m y of the apparent magni- 
tude of a quasar in two photometric bands x and y, we call 
C(y — x) = niy - m x the quasar color. If dust is present in a 
foreground absorption system, the observed color will be dif- 
ferent from the intrinsic, unreddened color Co(y - x), yielding 
a color excess E(y — x) — C(y — x) — Co(y — x). This cannot 
be measured directly since the unreddened color of the indi- 
vidual quasar is unknown. We can however measure the mean 
color excess of a sample of DLA-quasars in the following way. 
For each DLA-quasar we build a set of n i unreddened control 
quasars with same redshift and magnitude. We then estimate 
the mean color of these unreddened quasars, (C u (y - x)), and 
hence the deviation from the mean color 

A(y-x~) = C(y-x)- (C u (y-x)). (1) 

For a large sample of n, DLA-quasars we measure the mean 
color excess with the expression 

(E(y-x)) = (A(y-x)), (2) 



3.2. Choice of the color index 

Since our ideal goal is to measure the color of the 
quasar continuum, we excluded photometric bands contam- 
inated by the absorptions in the Ly a forest or the Ly 
a emission of the quasar. Taking into account the aver- 
age wavelengths A and the FWHMs of the SDSS photo- 
metric bands (Fukugita et al. 1996; see also SDSS web 
page http://www.sdss.org/dr5/instruments/imager) this condi- 
tion constrains the maximum quasar redshift to z fjma x = [(/l - 
iFWHM)/1216] - 1 = 1.67, 2.29, 3.50, 4.52, and 5.78 for the 
u, g, r, i, and z bands, respectively. Since the wavelength cov- 
erage of SDSS spectra prevents detection of DLA systems at 
z a < 2.2, only the r, i and z bands provide quasars with uncon- 
taminated photometry (in the sense specified above) and suf- 
ficient redshift coverage for detection of DLA systems. The 
maximum leverage for reddening detection with these bands 
is given by the (r - z) color. 

3.2.1 . The (r - z) color index 

In Fig.[T]we plot the color excess E(r-z) in the observer's frame 
that we expect to measure for an absorber at redshift z a with 
rest-frame extinction Ay - 1 magnitude. To estimate E(r - z) 
we considered different possible types of extinction curves for 
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Fig. 2. Mean color ((r - z)o) of the QSOs control samples used 
to estimate the unreddened colors of the DLA-QSOs. The mean 
color is plotted versus the QSO redshift, z e - Bright (r < 19.1 
mag) and faint (r > 19.1 mag) QSOs are plotted with blue 
crosses and red circles, respectively. 

the dust in the DLA systems, such as the Milky-Way curves by 
Cardelli et al. (1998), with the characteristic extinction bump 
at 2175A, and the SMC curves (Pei 1992, Gordon et al. 2003), 
characterized by a more regular UV rise. The extinction curves 
are normalized in the V band, %(X) = A A /A V , all the quantities 
in this definition being in the rest frame of the absorber. For an 
absorber at redshift z a the extinction measured at wavelength 

A° in the observer's frame is A^obs = Av^(jf^). The color 
excess in the observer's frame per unit extinction in the rest 
frame is therefore 



E(A° bs - X**) 

y — — =61;{z a ;Af\Af s ) 



(3) 



where 5%(z a \ /i? s , <tf s ) = - ^(jfr) ■ From thi s expres- 

sion we have computed the curves shown in Fig.Q] 

For an SMC-type extinction curve we expect a smooth vari- 
ation of the predicted color excess in the redshift interval of in- 
terest, with a typical value E(r - z) - 0.15 mag for absorbers 
with rest-frame Ay = 0.1 mag. For a MW-type curve we expect 
a small (or even negative) reddening in the interval of absorp- 
tion redshift of our sample. The minimum in the MW lines of 
Fig.fflcorresponds to minimum between the end of the extinc- 
tion bump at 2 175 A and the UV rise. The implications of these 
different behaviours of the SMC-type and MW-type extinction 
curves are discussed in Section 5.1. 



3.2.2. The (g - z) color index 

The measurement of the (g - z) color is challenging because 
the g band is generally contaminated by the Ly a forest, or 



even the QSO Lyman break, in the interval of emission redshift 
considered. However, measuring the reddening in the (g - z) 
color, in addition of the (r - z) color, offers several advantages. 
One is the higher leverage for reddening detection: we expect 
a gain of a factor of 2.2 in the ratio E(g - z)/E(r - z) at the 
redshift of our sample if the extinction curve is of SMC type. 
Another advantage is the possibility of probing the wavelength 
dependence of the reddening, i.e. the extinction curve of the 
dust, by comparing the (g - z) and (r - z) color excesses. We 
decided therefore to perform the measurements also in the (g - 
z) color index taking care of the contaminations present in the 
g bandpass. 

To avoid overlap of the g band with the QSO Lyman break, 
we restricted the analysis of the (g - z) colors to the QSOs 
with z e < 3.4 in our sample. We then took care of the contam- 
ination due to Ly a absorptions. The contamination due to the 
forest is present both in the quasars of the absorption sample 
and of the control sample. Line-to-line variations of the for- 
est absorption will contribute to the dispersion of the quasar 
colors, making the measurement of the color excess more dif- 
ficult. If such variations are randomly distributed, we still can 
measure the difference between the quasar colors of the ab- 
sorption and control samples. In doing this, however, we must 
take into account the fact that the g band of the absorption 
sample is generally contaminated by the damped Ly a profile, 
which is absent in the control quasars. To cope with this prob- 
lem we correct the g magnitude by taking into account the 
effect of the damped Ly a absorptiorQ in the filter bandbass. 
The magnitudes in the SDSS photometric system are defined 
as in — -2.5 log J <i(log v)f v S v + constant (Fukugita et al. 
1996), where f v is the quasar spectral distribution and S v the 
response function of the filter. We define the correction for the 
DLA absorption in the g band as 



5m„ 



-2.5 log 



fd(logv)f v °S 
jd(logv)fDS 



(a) 



(4) 



where f° is the quasar spectral distribution without DLA ab- 
sorption, ff the spectral distribution with DLA absorption, and 
S y the response function of the g filter (Fukugita, 2006, priv. 
comm.). To estimate 6m g we model f° with a power law with 
spectral index a v - -0.5, based on measurements of SDSS 
spectra of quasars in the redshift interval of interest for the 
present work (Vanden Berk et al. 2001, Desjacques et al. 2007). 
For each DLA system of our sample we then compute the the- 
oretical absorption profile of the lines of the Lyman series, t v , 
with intensity and position of the lines specified by the H i col- 
umn density and redshift. From this we compute /,f = l v X f° 
and hence the correction Sm g , which is finally added to the ob- 
served magnitude. The theoretical Voigt profiles are computed 
with the FITLYMAN routines (Fontana & Ballester 1995) in- 
cluded in the ESO MIDAS software package. The typical value 
of 5m g is of -0.05 mag, with values as low as -0.006 mag 
and as high as -0.17 mag in the most extreme cases. 



4 A similar procedure was adopted by Ellison et al. (2005) to cor- 
rect for the flux suppression due to Ly a absorption in the B band. Here 
we consider the flux suppression due to all the lines of the Lyman se- 
ries that overlap the SDSS g band. 



Vladilo, Prochaska, & Wolfe: color excess of QSOs with intervening DLA systems 



7 



3.3. Implementation of the method 

The "BEST" ugriz imaging data were recovered for all quasars 
of the absorption and control samples. These photometric data 
were corrected for the effects of the Galactic extinction accord- 
ing to the prescriptions given by Schneider et al. (2005). 

For each quasar of the absorption sample we calculate the 
colors (r - z)/ and (g - z)i and then select a taylored subset 
of control quasars closest in redshift and magnitude in the fol- 
lowing way. We first select the n Ze control quasars closest to 
the redshift z e of the DLA quasar. From this subset we then 
select the n m control quasars closest to the z-band magnitude 
m z (the least affected by extinction) of the DLA quasar. This 
last subset is used to estimate the mean unreddened colors 
(i r - Z)o)i and (ig - z)o);- The deviations from the mean color 
are finally calculated as A,(r — z) = (r - z)i - ((r - z)o)i and 
A,-(g -z) = (g- z)i - ((g - z)o)j. 

With this procedure each subset of control quasars has the 
same size n m for all DLA-QSOs and we have a comparable 
statistics for all the measurements. The choice of n Ze and n m is 
determined by the requirement that the total intervals in emis- 
sion redshift, Si(z e ), and z magnitude, Si(m z ), spanned by each 
subset are sufficiently small to guarantee a good degree of ho- 
mogeneity of the control samples. In practice, this gives the 
contraints n Ze < 100 and n z Jn m > 4 in order to have typical 
(median) values (5,(z e ) ^0.1 and <J,(m z ) 0.5 mag. The red- 
shift interval is similar to the redshift bins commonly adopted 
in studies of the mean quasar color as a function of z e (e.g. 
Richards et al. 2001). Changes of the intrinsic slopes of the 
quasar continua over the above magnitude interval are expected 
to be modest (Yip et al. 2004). 

In Table 1 we report the measurements of A,(r - z) and 
- z) for the DLA-QSOs of the absorption sample. In 
columns 6 and 7 of the table we give the intervals in emis- 
sion redshift, Si(z e ), and z magnitude, <5,(m z ), spanned by each 
subset of control quasars obtained by adopting n Zt = 100 and 
n m = 25. In columns 8 and 9 we list the mean colors {{r - z)o)i 
and ((g - z)o)j, obtained from the weighted mean of the colors 
of each control subset. The weights are taken to be proportional 
to the inverse squares of the errors of individual colors. 

The errors of individual colors of DLA quasars and control 
quasars are obtained by propagation of the photometric errors 
and of the errors of the Galactic extinction correction. For the 
latter we adopted an uncertainty of +50% of the correction it- 
self. The (g - z) data include the correction for damped Ly a 
absorption in the g band, when approprite. An uncertainty of 
±50% of this correction is propagated in the error budget. 

The errors of the mean colors {{r - z)o>; and ((g - z)n); 
quoted in Table 1 are the standard deviation of the colors of 
each subset. These errors dominate the budget of the errors of 
the color deviations A,(r - z) and A,(g - z). 

In Fig.|2]we plot the mean colors ((r - z)o), of the subsets of 
unreddened quasars versus z e . Each data point corresponds to 
the subset of an individual DLA-quasar of Table 1 . We use dif- 
ferent symbols for quasars brighter and fainter than the median 
magnitude of the absorption sample, m r = 19.1. In addition to 
a trend with z e , known from previous studies, the figure shows 
a dependence on the quasar magnitude. For instance, at z e ~ 3, 



where many quasars both fainter and brighter than m r = 19.1 
are present, the fainter quasars (red circles) lie systematically 
below the brighter ones (blue crosses). This result justifies our 
choice of grouping the control quasars not only in redshift, but 
also in magnitude. 

4. The mean color excess 

In Table 3 we give the weighted means (E(r - z)) and 
(E(g - z)> estimated with the expression (E) = (A) = 2/ vv/A,-, 
where w, = o"7 2 / 2, o~ 2 , and <x, are the errors of individual 
measurements. The use of the weighted mean allows us to min- 
imize the contribution of the QSOs with largest uncertainty of 
their intrinsic colors, the dominant source of the error budget. 
For the error of the mean we adopted the unbiased estimate 



2 / w / (A i -<A»2 



(n-l) 



(5) 



(Linnik 1961) with the normalized weights w, given above. 
This expression is valid if the measurements of A, are uncor- 
rected. We performed two tests to make sure that correlated 
errors are unimportant and therefore the adopted statistical er- 
ror cr m is reliable. First we computed the covariance matrix and 
compared the strength of off-diagonal terms with that of diago- 
nal terms. The covariance matrices of the A,(r - z) and A/(g - z) 
data samples were computed using a bootstrap method over 60 
random realizations of each sample. We found that the normal- 
ized covariance matrices are fairly close to be diagonal, with 
off-diagonal terms tipically around a few percent, and always 
below 25%. The low values of the off-diagonal terms indicate 
that there is little correlation between the measurements. Then 
we studied the behaviour of <x„, as a function of the sample 
size. This can be done with a random grouping technique de- 
scribed by Hampel et al. (1986; Chap. 8). From this test we find 
that the adopted error follows the theoretical n~ decrease, as 
expected for uncorrected measurements. The adopted error of 
the mean is therefore a good estimate of the statistical error of 
(E). 

4.1. Confidence level of detection 

From the results shown in Table 3 one can see that the mean 
(r - z) reddening is detected at ^ 3 cr m confidence level (statis- 
tical error) in the complete sample of 248 DLA/QSO pairs for 
which the r band is uncontaminated by Ly a forest (z e < 3.5). 
The mean (g - z) reddening is detected at ^ 4.5 cr m confidence 
level in the sample of 232 DLA/QSO pairs with g magnitude 
uncontaminated by the QSO Lyman break (z e < 3.4) and cor- 
rected for DLA flux suppression (Section 3.2.2). 

To estimate the confidence level of the detection other than 
in terms of cr m we applied the method of bootstrap resampling 
(Efron 1979). A bootstrap sample is obtained by extracting n 
data at random with repetition from the original sample of n 
measurements. This type of extraction is, in practice, an ex- 
traction with replacement because part of the original measure- 
ments is replaced by repeated data. By iterating this extrac- 
tion process many times it is possible to build a large number 
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Table 3. Mean reddening in the observer's frame, (E), and mean extinction in the rest frame, (Ay), for different colors and 
DLA/QSO samples. 
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Nb of bootstrap samples from the original set of n data. For 
each of the Nb samples one can compute a statistical estimator 
of interest (e.g. the weighted mean) and study the distribution 
of such estimator, without making any assumption on the par- 
ent distribution. We applied this method to the original sample 
of n measurements A,, from which we built up Nb = 10,000 
bootstrap samples. The weighted mean (A) was recomputed in 
each case. We analyzed the resulting distribution of (A) val- 
ues (Fig.0 to estimate the fraction of cases in which the mean 
value is positive. We find (A(r - z)) > in a 99.90% of the 
_ z) bootstrap samples and (A(g - z)) > in a 99.99% of 
the Aj(g - z) bootstrap samples. These percent figures indicate 
the probability that a positive color excess has been detected. 

4.2. Analysis of "twin control samples" 

A reason of concern in the present analysis is the possibility 
that the mean color excess that we detect is due to sources of 
reddening other than the DLA systems (e.g. dust in the quasar 
environment or in unidentified absorbers). If present, these red- 
dening contribution would also affect the quasars of the control 
sample. To understand if these effects may accidentally yield 
a mean color excess comparable to that measured in the DLA 
sample, we applied our procedure to a large number of "twin 
control samples" with same size, redshift distribution and mag- 
nitude distribution as the DLA sample. We then analyzed the 
values of (A) obtained from these twin control samples in a 
search for cases with a value equal to those reported in Table 3. 
In practice, each twin control sample was created by extracting 
at random, for each DLA-QSO of Table 1 , a control QSO at the 
same redshift and magnitude. The control QSOs extracted in 
this way were treated as DLA-QSOs and temporarily excluded 
from the estimate of the mean unreddened colors (C) (Section 
3.1). As a consequence, the estimates of (C) were different for 
each realization of the twin control sample. We built up 10,000 
twin control samples and computed the weighted mean (A) for 
each of them. The resulting frequency distribution of (A) val- 
ues is shown in Fig. [4] for both color indices of interest. For 
comparison we show the mean color excess of the absorption 
sample (dotted lines). The probability for a twin control sam- 



ple to attain a value (A) as high as the mean color excess of the 
absorption sample is < 0.01% in (g - z) and a 0.4% in (r - z). 
We conclude that the mean color excess detected in the absorp- 
tion sample cannot be ascribed to dust in the quasar environ- 
ment or in unidentified low-redshift absorbers. The same low 
probability applies to any sort of systematic effect that might 
accidentally conspire to yield a mean color excess as a conse- 
quence of the application of our procedure to a quasar sample 
with same size, redshift distribution and magnitude distribution 
as the DLA sample. 



4.3. Rejection of the most reddened quasars 

For a correct interpretation of the present results it is important 
to understand whether the mean reddening that we detect is a 
signature of the bulk of the DLA population or is just due to 
the presence of a few reddened quasars. To clarify this point 
we rejected the most reddened DLA-QSOs from the original 
sample of Table 1 and repeated the computation of the mean 
reddening. In practice, we discarded the cases with individual 
color deviation A,- larger than the mean color deviation (A) at 
5 cr level ("absorption sample 2") and at 3 cr level ("absorption 
sample 3"). In the second and third row of Table 3 we show 
the mean color excess recomputed after rejecting the most red- 
dened cases in this way. One can see that the mean color excess 
is still detected at a 2.8 cr m level in (r - z) and at 3.5 cr m level 
in (g - z) for the absorption sample 3. The detection is firmer 
for the absorption sample 2. 



4.4. Size of the subsets of control quasars 

As an final test we checked the stability of the results for the 
adoption of different values of n Ze and n m (Section 3.3). We 
tested this effect by varying these parameters around their opti- 
mal values n Ze a 100 and n z Jn m a 4. The resulting changes of 
the mean reddening are significantly smaller than the statistical 
error discussed above. 
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Fig. 3. Frequency distribution of the mean color excess (A(r - z)) and (A(g - z)) measured in 10,000 bootstrap samples obtained 
from the list of A values of Table 1 . Dotted line: mean color excess of the original sample. See Section 4. 1 . 
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Fig. 4. Frequency distribution of the mean deviations (A(r - z)} and (A(g - z)} measured in 10,000 control samples with identical 
size, redshift distribution and magnitude distribution of the absorption sample. Dotted line: mean color excess of the absorption 
sample. See Section 4.1. 



5. Discussion 

The analysis presented in the previous section indicates that 
high-redshift quasars with intervening DLA systems at 2.2 < 
z a < 3.5 have redder colors than quasars at the same redshift 
without foreground absorption systems. 

The most natural process to explain the origin of this red- 
dening is dust extinction, with its characteristic rise at shorter 
wavelengths. If this is the case, we expect the measured red- 



dening to show similarities with that predicted for typical in- 
terstellar dust extinction curves. We also expect the most red- 
dened quasars to be dimmed by dust extinction and therefore 
to be statistically fainter than unreddened quasars. In Sections 
5.1 and 5.2 below we present observational evidence consistent 
with these expectations. 

The location of the dust is of special interest in the present 
work. Reddening sources in the quasar environment cannot ex- 
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plain the difference in color between the quasars of the absorp- 
tion and control samples since this type of reddening would af- 
fect in the same way, on the average, the quasars of both sam- 
ples. Dust embedded in the DLA systems is the most natural 
hypothesis to explain the measured color excess. This hypoth- 
esis bears some predictions that can be tested empirically. 

One is that the reddening should increase with the column 
density of the metals embedded in the DLA systems. Another is 
that the dust-to-gas ratio should approximately scale with the 
level of metallicity of DLA systems. In Sections 5.3 and 5.4 
below we present observational evidence consistent with these 
expectations. In Sections 5.5 and 5.6 we discuss the reddening 
versus H i column density and versus absorption redshift. 

5.1. Dust extinction curve 

One of the properties used to characterize the interstellar dust 
is its extinction curve, g(A) = A^/Ay. In spite of strong spatial 
variability, interstellar extinction curves can be classified in two 
types: Milky-Way type curves, with the characteristic extinc- 
tion bump at 2175 A (e.g. Cardelli et al. 1998), and SMC-type 
curves, characterized by a more regular UV rise (e.g. Pei 1992, 
Gordon et al. 2003). The different slopes of these curves can be 
used to distinguish between them. 

If the reddening originates in the DLA systems, we can 
measure the slope of the DLA dust extinction curve by com- 
paring the mean reddening in the (g - z) and (r - z) colors. We 
use the restricted sample of 232 DLA-QSOs with z em < 3.4, 
for which both (g - z) and (r - z) can be derived. For this 
sub-sample we measure (E(r - z)) = 25(±9.0) 10~ 3 mag and 
we obtain (E(g - z))/(E(r - z)) = 2.2 + 0.9. This measure- 
ment is in good agreement with the ratio predicted for ab- 
sorbers with SMC-type extinction curve at 2.2 < z a < 3.5, 
(E(g - z) I 'E{r-z)f^ -2.2. 

The fact that we detect reddening in the (r - z) color argues 
againts a MW-type extinction curve. In fact, this type of curve 
is expected to give a negligible or even negative reddening in 
the interval of absorption redshift spanned by of our sample 
(Fig. [TJ. This is due to the fact that in this interval the (r - z) 
color probes a part of the MW extinction curve with negative 
slope, namely the end of the 2175 A extinction bump. The fact 
that we detect reddening in the (r - z) color indicates that ex- 
tinction curve of MW type must be uncommon in the sample. 

If DLAs with MW-type curves exist, a broad absorption 
centered at (1 +z fl )x2175 A should affect the quasar continuum. 
We searched for this feature in all the spectra of the absorption 
sample, and in particular of the quasars with anomalous colors, 
but without success (see e.g. Fig. [SJ- 

The above results are in line with the general finding that 
MW-type extinction curves are rare among quasar absorbers 
(Junkkarinen et al. 2004, Wang et al. 2004), while SMC-type 
are common (Wild & Hewett 2005, York et al. 2006). By adopt- 
ing an SMC extinction curve, our results imply a rest-frame ex- 
tinction (Ay) - 0.015/0.019 mag (Table 3), corresponding to 
(E(B - V)> - 0.005/0.007 mag. 



5.2. Reddening versus quasar magnitude 

As we mentioned above, we expect the most reddened quasars 
to be statistically fainter than unreddened quasars if the red- 
dening originates in dust. In Fig. [5] we plot A(r - z) and 
A(g— z) versus the r magnitude. The largests values of color ex- 
cess lie at faint magnitudes (r > 18 mag), in broad agreement 
with this expectation. To quantify the effect we compared the 
mean magnitude of the 30 DLA-QSOS with highest A(r - z) 
with the mean magnitude of the remaining DLA-QSOs with 
lower A(r - z). The most reddened cases are slightly fainter 
on the average, (r) - 19.24 ± 0.08, than the remaining cases, 
(r) = 18.98 + 0.04. A similar computation for the (g - z) color 
yields (r) = 19.11 ± 0.09 for the 30 most reddened cases and 
(r) = 18.96 + 0.05 for the remaining cases. The effect is mod- 
est, but is in line with the prediction that quasars with higher 
reddening should be statistically fainter. 



Table 4. Mean reddening of sub-samples with different 
strengths of the Sin 1526 line detected at the redshift of the 
DLA system. 







VVsin 


< 1 A 




> 1 A 


color 


Sample 


n 


(E) 


n 


(E) 








[10~ 3 mag] 




[10- 3 mag] 


(r-z) 


1 


151 


26 ± 12 


29 


41 ±25 




2 


150 


21 ± 11 


29 


41 ±25 




3 


149 


17 ± 10 


29 


41 ±25 


(g-z) 


1 


142 


51 ± 14 


29 


114 ±40 




2 


141 


46 ± 13 


29 


114 ±40 




3 


140 


42 ± 13 


29 


114 ±40 



5.3. Reddening versus metal column density 

If the reddening originates in dust embedded in the DLA sys- 
tem we expect a trend between the color excess and the total 
column density of metals in the absorber. The limited spec- 
tral resolution of SDSS spectra prevents the accurate mea- 
surement of metal column densities, but is sufficient to mea- 
sure the equivalent width of strong metal lines. To search for 
a possible trend between color excess and metals we per- 
formed a study of the Sin line at 1526 A, which is suffi- 
ciently strong to be detected in SDSS spectra and lies red- 
wards of the Ly a forest in most of the redshift range of in- 
terest. Since the 1526 A line is saturated, its rest-frame equiva- 
lent width, Wsi n, is not a good tracer of the metal column den- 
sity. Nevertheless, an empirical correlation between Wsin and 
metallicity has been recently found for DLA-QSOs, in the form 
[M/H] - -0.92 + 1.41 log Wsi„ (Prochaska et al. 2007). This 
correlation is probably the consequence of the mass-metallicity 
and mass-kinematics relations in galaxies. Galaxies with higher 
masses (and metallicities) have larger kinematical motions, 
which result in larger spread (and equivalent width) of satu- 
rated metal lines such as the Si n line at 1526 A. 

To estimate Wsin we used an automated algorithm which 
fits a Gaussian profile to the strongest absorption feature (if 
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Fig. 5. Deviation from the mean color, A(r - z) and A(g - z), versus r magnitude for sample 2 of Table 3. Dotted lines: mean 
color excess of the same sample. 
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Fig. 6. Deviations from the mean color, A(r - z) and A(g - z), versus rest-frame equivalent width of Sin line 1526 A, when 
detected, for sample 2 of Table 3. Dotted lines: mean color excess of the same sample. 



any) at the predicted wavelength, Asin = (1 + z a ) 1526.707 A, 
provided this line occurs redwards of the Ly a forest. Each fit 
was visually inspected and the line-strength of Sin 1526 was 
compared against other metal-lines observed redwards of the 
Lya forest. In a few cases, we rejected lines because of ob- 
vious blends with coincident absorption lines from systems at 
unrelated redshifts. As a result of this analysis, we measured 
an equivalent width in 180 DLA-QSOs of our sample. In Fig. 
[6] we plot the deviations from the mean colors A(r - z) and 



A(g - z) versus the rest-frame equivalent width Wsin measured 
by the automated algorithm. No clear trend is present, possibly 
because most of the data are clumped at low values of equiva- 
lent width, where we expect a modest reddening. The interval 
Wsin > 1 A, where the reddening is expected to be stronger, is 
much less populated, but its size is sufficiently large to compute 
the mean value. We therefore splitted the data in two subsets 
of "weak" and "strong" Sin absorbers by using Wsin = 1 A 
as a threshold value. The mean values of reddening derived for 
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these subsets are shown in Table 4 for the three absorption sam- 
ples defined in Table 3. One can see that the mean reddening 
tends to be higher in strong Si n absorbers than in weak Si n ab- 
sorbers, as expected if the dust is embedded in DLA systems. 
The difference in mean reddening between the two sub-samples 
is marginal in the (r - z) color, but more evident in the (g - z) 
color. 

To investigate the physical relation between dust and metals 
we converted the color deviations A into a mean extinction and 
the Si ii equivalent widths into a metallicity. The rest-frame ex- 
tinction, Ay, is expected to scale with the metal column density 
(M/H) x N(H i), where (M/H) is the total abundance by number 
(gas plus dust) of the reference element M (e.g. Vladilo et al. 
2006). 

To derive the extinction we used the expression Ay = 
E(A° bs - Af s )/6^(ZaUf s ,A° hs ), obtained from Eq. ©. The 
true color excess Ej of the individual quasar is unknown, but 
we can take E, A, since we average the above expression 
for a large sample to obtain (Ay) (see Section 3.1). An SMC 
extinction curve £(A) was adopted in this conversion. 

To derive the metal column density we used the empirical 
relation [M/H] -0.92 + 1.41 log W si „ (Prochaska et al. 2007) 
and then averaged the individual values of metal column den- 
sity lolWfU+ketfCHi^ 

In Fig.|7]we plot the results of these computations obtained 
for the subsets of weak and strong Si n absorbers of Table 4 
(sample 2). The data are consistent with the existence of the 
trend between metal column density and extinction expected if 
the reddening and metals originate in the same medium. Again, 
the evidence for this trend is very weak in (r - z) and firmer in 
(g-z). Given the uncertainties in the above conversions and the 
limited statistics it is not possible to reach firmer conclusions. 

The quasars of the absorption sample with high color devi- 
ations and strong Si II absorptions are excellent candidates for 
detailed follow up studies of dusty DLA systems based on high 
resolution spectroscopy. In Fig.|8]we show, as an example, the 
spectra of the DLA-quasars with A{g-z) > 0.4 and Wsin > 1 A. 
The typical metallicity of these systems is a factor 3 larger than 
the average metallicity of the absorption sample. 

5.4. Dust-to-gas ratio 

By adopting an extinction curve £(A) of SMC-type we can con- 
vert the reddening measured in the observer's frame into an 
extinction in the rest frame. From this in turn we can measure, 
for the first time in studies of DLA systems, the dust-to-gas ra- 
tio in the same units used in local interstellar studies, namely 
the ratio Av/N{Hi) between the extinction in the V band and 
the Hi column density. From Eq.® the individual dust-to-gas 
ratio is given by 



AT(Hi) 



E(A™-A™)/6Z(zgUT>AT) 
JV(Hi) 



(6) 



Since we average this quantity for all the systems of our sam- 
ple we take Ei - A, (see Section 3.1). The mean values ob- 
tained from this computation are shown in Table 3. We obtain 
(Av/N(Hi)) ^ 3 to 4 xl0~ 23 mag cm 2 , from the measurements 
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Fig. 7. Mean values of rest-frame extinction and metal col- 
umn density of the two subsets of DLA-QSOs with different 
strengths of the Sin line 1526 A shown in Table 4 (absorption 
sample 2). The extinction values derived from the (g - z) and 
(r - z) color deviations are plotted with diamonds and circles, 
respectively. The mean values of metal column density of the 
two subsets are the same and have been slightly shifted along 
the horizontal direction for clarity. See Section 5.3. 



in the (r—z) or (g-z) colors, respectively. By excluding the most 
reddened cases we obtain (Ay/N(Hi)) 2.5 to 3.5 10 -23 mag 
cm 2 . For comparison, the typical dust-to-gas ratio of the Milky 
Way is A V /N(H i) * 5.3 x 10~ 22 mag cm 2 (Bohlin et al. 1978). 
These results indicate that the dust-to-gas ratio of DLA systems 
is deficient by » - 1 .25 dex relative to that of the Milky Way. 
This deficiency is consistent with the lower level of metallicity 
of the DLA systems of our sample. This metallicity can be es- 
timated from the relation between [M/H] and Wsin (Prochaska 
et al. 2007). For the systems with measurements of Wsin we 
obtain a mean value ([M/H]) - —1.2 dex (Table 5). We con- 
clude that the dust-to-gas ratio scales with the metallicity, as 
expected for dust embedded in the DLA galaxies. 

5.5. Dust-to-metal ratio 

The present measurements of reddening and metallicity can be 
combined to estimate the mean extinction per unit metal col- 
umn density in the DLA systems of our sample. We start from 
the relation 



<4?>x/ Fe x(Fe/H)xjV(Hi) 



(7) 



where (Sy e ) is the mean extinction per unit dust column density 
of iron and fy e is the fraction of iron in dust form (Vladilo et al. 
2006). Previous measurements of (sy e ) have yielded a typical 
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Fig. 8. Quasars of the absorption sample with color deviations A(g — z) > 0.4 mag and strong Si n absorptions at zsm = Zdla- 
Dashed line: profile of the 2175 A extinction bump plotted in a scale l(H [i/(1+ZDLA)] (see Vladilo et al. 2006). 



Table 5. Mean metallicity and extinction per unit metal column density for the DLA-QSOs with Si n measurements. 



color 


(r-z) 


(g-z) 


Absorption 


n 


<[M/H]) 6 






n 


<[M/H]> 6 


<[m/h]>;;. 




sample" 








[10~ 17 mag cm 2 ] 








[10~ 17 mag cm 2 ] 


1 


180 


-1.19 ±0.03 


-1.07 


14 ±8 


171 


-1.18 ±0.03 


-1.07 


7.5 ± 4.3 


2 


179 


-1.18 ±0.03 


-1.07 


12 ±8 


170 


-1.18 ±0.03 


-1.06 


5.8 ±4.0 


3 


178 


-1.18 ±0.03 


-1.07 


11 ±8 


169 


-1.17 ±0.03 


-1.06 


5.4 ± 4.0 



" Absorption samples are defined in Table 3. Only cases with Sin line detected at more than 1 cr level are considered in this table. 
* Mean metallicity {[M/H]> = log (10 [M ' H) ). The quoted error is the standard error of the mean. 
c Mean Hi column-density weighted metallicity of Eq. 

d Mean extinction per unit dust column density of iron (Section 5.5). The quoted error is the standard error of the mean. 



valu^fl (Sy e ) « 3 x 10~ 17 mag cm 2 , roughly constant within 
a factor of two, in different types of Hi interstellar regions, 
including metal absorption systems and a few DLA systems at 
z a < 2. The present data set allows us to investigate whether 
or not this parameter has a similar value also in the general 
population of DLA systems at z a > 2. 

By combining the Eqs. (0 and (0 we obtain 



<4 e ) = — 



A? s )/5t(z a ;<i° 



/ Fe x(Fe/H)xjV(Hi) 



(8) 



As in the previous section, we take E-, - A,- because we av- 
erage this expression for all the systems of our sample. We 
adopt an SMC extinction curve, assume (Fe/H) = (M/H) and 
take /p e - 0.65, a value appropriate for DLA systems at the 
redshift of our sample (Vladilo 2004). The resulting values of 
(iy e ), shown in Table 5, are broadly consistent with the value 
(iy e ) » 3x 10~ 17 mag cm 2 inferred from previous investiga- 
tion^] at lower redshift. This is true, in particular, for our best 



5 Mean value of the most accurate measurements of A v and dust 
column density Nf c for the DLA systems in Table 3 by Vladilo et al. 
(2006); the same value matches well the interstellar trend between N Fs: 
and A v (see Fig. 4 in the same paper). 

6 The recent result (E(B - V))/(N Zn ) ^ 4.1 x 10~ 15 mag cm 2 , found 
by Menard et al. (2007) in Mgn systems, yields the same value of 

assuming that the Zn/Fe ratio is solar and /p e ^ 0.65. 



estimate obtained from the (g-z) color deviations. Considering 
the uncertainties in the values of fy e and (Fe/H), this broad 
agreement is rather remarkable and suggests that the extinction 
properties of the dust are similar to those of low -redshift dust. 
A firmer conclusion must await more accurate determinations 
of the metallicities and depletions of the present sample. 



5.6. Reddening and H i column density 

In Fig.|9]we plot the color deviations A(r-z) and A{g-z) versus 
H i column density. A linear correlation analysis of A(r - z) and 
A{g - z) versus log A^(H i) for the complete absorption sample 
yields correlation coefficients r = -0.01 and r = +0.01, re- 
spectively. The comparison of sub-samples with low and high 
values of A^(H i) does not indicate the existence of a trend with 
reddening. Most DLA-QSOs are concentrated at low Hi col- 
umn densities as a consequence of the well-known decrease of 
the number of DLA systems with increasing A^(H i) (Wolfe et 
al. 2005; Prochaska et al. 2005). Some of the most reddened 
quasars are concentrated at low Hi column densities, at odd 
with the expectation that large reddening is associated with 
large column density of neutral gas. The paucity of highly red- 
dened cases at A^(Hi) > 10 21 atoms cm~ 2 could be the con- 
sequence of low-number statistics: systems with high Hi col- 
umn density are rare and the present statistics might be insuffi- 
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Fig. 9. Deviations from the mean color, A(r - z) and A(g - z), versus Hi column density for sample 2 of Table 3. Dotted lines: 
mean color excess of the same sample. 



cient to detect cases of high reddening. It is worth noticing that 
lines of sight of high metallicity are generally not observed in 
DLA systems with N(Hi) > 10 21 atoms cnT 2 (Boisse et al. 
1998). If this is the case also for the present sample, this fact 
could conspire to make rare the systems of high reddening at 
N(Ui) > 10 21 atoms cirT 2 . 

5. 7. The H i column-density weighted metallicity 

The metal abundances of DLA systems are often used to es- 
timate the mean cosmic metallicity of the neutral gas at high 
redshift. In this type of estimate the metallicity of each system 
must be weighted by its H i column density to estimate a mean 
cosmic value 

2,(M/H),xM(Hi) 



<[M/H]> B . 



(9) 



(Lanzetta et al. 1995). We estimated the above expression by 
using the indirect measurement of metallicity obtained from 
the equivalent width of the Sin line. The results are shown in 
Table 5. We obtain ([M/H]) w - -1.1 dex. This value is some- 
what higher than the mean Hi- weighted metallicity of DLA 
systems measured in high-resolution surveys at redshift z a ~ 3, 
which is s —1.5 dex (Prochaska et al. 2001; Kulkarni et al. 
2005). The present survey of low -resolution SDSS quasar spec- 
tra is characterized by a higher magnitude limit since it includes 
quasars as faint as m r = 20.2. The higher value of metallicity 
of the present SDSS sample suggests that the mean metallic- 
ity of DLAs may increase with the limiting magnitude of the 
survey. An effect of this type is predicted as a consequence of 
the dust extinction bias (Ellison et al. 2001, Vladilo & Peroux 
2005), but the empirical evidence for its existence is still lack- 
ing (Akerman et al. 2005). An accurate determination of the 
metallicities of the present absorption sample will offer us a 



way to test the dependence of the mean metallicity on the lim- 
iting magnitude of DLA-QSO surveys. 



5.8. Reddening versus DLA redshift 



In Fig. [TO] we plot the color deviations A(r - z) and A(g - z) 
versus absorption redshift z a . The data are spread more or less 
uniformly in the interval of absorption redshift considered. A 
linear correlation analysis of A(r - z) and A(g - z) versus z a for 
the complete absorption sample yields correlation coefficients 
r - +0.13 and r = +0.06, respectively. The lack of a trend 
with redshift implies that the color excess that we derive is suf- 
ficiently representative of the mean reddening of DLA systems 
over the full interval 2.2 < z a < 3 .5 covered by our sample. The 
lack of redshift evolution is, to some extent, surprising. Two 
types of evolutionary effects are expected to occur if the red- 
dening originates in DLA dust. On the one hand, a rise of the 
reddening with increasing z a is expected due the UV rise of 
the extinction curve. For an SMC extinction curve this effect 
is roughly a factor of 2 in E(r - z) in the redshift range of in- 
terest, as shown in Fig. [JJ On the other hand, a decrease of the 
reddening with increasing z a is expected as a consequence of 
cosmic chemical evolution if the dust scales with the metals. 
Also this effect is expected to be approximately a factor of 2 
in the redshift interval of interest, where the mean metallicity 
decreases by 0.3 dex per unit redshift interval (Prochaska et 
al. 2000). The lack of a trend could be the consequence of the 
opposite behaviour of two evolutionary effects, which tend to 
cancel each other. 
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Fig. 10. Deviations from the mean color, A(r - z) and A(g - z), versus absorption redshift for sample 2 of Table 3. Dotted lines: 
mean color excess of the same sample. 



5.9. Comparison with previous work 

To perform a homogeneous comparison with the work of 
Murphy & Liske (2004; ML04) we applied our procedure to 
their original list of quasar^] extracted from the second data re- 
lease (DR2). The list includes 70 DLA-quasars and 1396 non- 
DLA quasars in spectra with SNR > 3 as defined in ML04. 
To apply our procedure we selected the 56 DLA-quasars and 
1074 non-DLA quasars with z e < 3.5. The non-DLA quasars 
were used to build the control pool. From the application of 
our method to these lists we obtain the weighted mean val- 
ues (E(r — z)) = 3 (+22) x 10~ 3 mag and, for the subset of 
52 quasars with z e < 3.4, {E(g-z)) = 39 (+36) x 10~ 3 mag. 
The lack of detection of the reddening is due to the small size of 
the sample. By converting these results to rest-frame color ex- 
cess with an SMC extinction curve we obtain a 3<x upper limit 
E(B - V) < 0.02 mag from both color indices. This limit is 
equal to that obtained by ML04 from their analysis of the spec- 
tral index distribution of the quasar spectra. The photometric 
method presented here and the spectral index method adopted 
by ML04 yield therefore consistent results. 

6. Summary and conclusions 

We have used the spectroscopic and photometric database of 
the 5th data release of the SDSS to measure the reddening of 
quasars with intervening Damped Ly a systems at high red- 
shift (zdla > 2.2). The spectroscopic database was used to dis- 
tinguish quasars with and without intervening absorption sys- 
tems. Only good quality spectra, with typical signal-to-noise 
ratio SNR > 4 in the spectral regions of interest, were con- 
sidered in the analysis. We built up an "absorption sample" of 

7 See web address http://www.blackwellpublishing.com/products 
/journals/suppmat/MNR/MNR8374/mnr8374TableS 1 .txt. 



248 quasars with a single DLA system in the redshift interval 
2.2 < z a < 3.5 and without metal systems at z me tais + Zdla- 
We then selected about 2 thousand control quasars without 
DLA systems or metal systems. The large SDSS photomet- 
ric database, employed for the first time in the DLA reddening 
analysis, was used to measure the colors of individual quasars 
of the absorption and control samples. 

The analysis was performed on the (r - z) and (g - z) color 
indices. The color index (r - z) is uncontaminated by Ly a ab- 
sorption in the redshift interval considered. The color (g - z) 
offers a higher leverage for reddening detection, but the g band- 
pass is contaminated by Ly a absorption. To derive the {g - z) 
colors we corrected the g magnitude for flux suppresion due to 
the presence of the damped Ly a absorption in the spectra. We 
obtain a mean color excess (E(r-z)) = 27 (+9) X 10 -3 mag and 
(E(g - z)) = 54 (±12)xl0~ 3 mag. The quoted statistical erroris 
not significantly affected by the presence correlated errors. The 
detection is confirmed by the analysis of 10,000 bootstrap sam- 
ples originated from the absorption sample. The reddening is 
detected at » 3 <x c.l. also when the most reddened DLA-QSOs 
of the sample are rejected, assuming that they are spurious. 

The most natural hypothesis to explain the observed red- 
dening is dust along the line of sight. The ratio (E(g-z))/(E(r- 
z)) = 2.2 + 0.9 that we measure lends support to this hypoth- 
esis, since this value is in agreement with that predicted for 
dust with SMC-type extinction curve at the absorption redshift 
of our sample. Also the study of the reddening versus quasars 
magnitude is consistent with the dust hypothesis: the most red- 
dened quasars are, on the average, slightly fainter than the oth- 
ers, in line with the expectation that extinction and reddening 
should increase together if the reddening is due to dust. 

Dust in the quasar environment or in low-redshift interlop- 
ers is not a viable explanation for the observed reddening since 
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it would affect in the same way, on the average, the quasars of 
the absorption sample and those of the control sample. Quasars 
with anomalous colors would not be able to produce acciden- 
tally the measured reddening (Section 4. 1). 

The above arguments indicate that the reddening is due 
to dust embedded in the intervening DLA systems. This con- 
clusion is consistent with two other studies performed in the 
present work. One is the comparison of the color excess with 
the equivalent width of the Sin 1526 A line at the redshift of 
the DLA system: strong Sin absorbers (Wsin > 1 A) show a 
slightly higher reddening, on the average, than weak Si n ab- 
sorbers (Wsin < 1 A), in line with the expectation that the 
amount of dust and metals should increase together inside the 
DLA systems. Also the study of the mean dust-to-gas ratio, 
obtained by converting the observed reddening into rest-frame 
extinction, yields a similar conclusion. We find (A v /N(Hi)) « 
2 to 4 Xl0~ 23 mag cm 2 , a value * -1.25 lower than that of 
the Milky Way, consistent with the lower level of metallicity 
of DLA systems of our sample estimated indirectly from Wsin. 
Also this result is in line with the expectation that metals and 
dust should increase together. 

The conversion of our reddening measurement to rest- 
frame color excess by means of an SMC extinction curve yields 
(E(B - V)> - 0.005/0.007 mag. This value fits the trend of 
{E(B - V)) versus redshift in quasar absorbers extrapolated 
from measurements of Mg n systems at lower redshift (Menard 
et al. 2007; Fig. 10). 

The mean color excess that we have derived is represen- 
tative of the mean reddening of DLA systems at z a ~ 2.7 in 
SDSS QSOs with limiting magnitude r 20.2. The magnitude 
limit results from the condition SNR > 4, required to distin- 
guish the presence of absorption systems in the quasar spectra. 
Since highly reddened quasars are more likely to be found at 
fainter magnitudes, a larger reddening could be obtained from 
deeper surveys. 

Follow-up, direct measurements of metallicity of the DLA 
systems of the present sample are required to test the indirect 
estimates based on W$i U , to obtain an accurate measurement 
of the extinction per unit metal column density, which seems 
to be similar to that found at lower redshift, and to assess the 
importance of the extinction bias in studies of DLA-QSOs. 

Acknowledgements. We thank Gabriel Prochter for providing an up- 
dated list of Mgn absorptions in SDSS DR5 quasar spectra. G.V. 
thanks Matteo Viel, Sergei Levshakov and Carlo Morossi for help- 
ful discussions on the analysis of the data. J.X.P. is partially supported 
by an NSF CAREER grant (AST-0548180). We thank an anonymous 
referee for comments that improved this paper. 

References 

Adelman-McCarthy, J. K., Agiieros, M. A., Allam, S. S., Kurt S. J., 
Anderson, K. S. J. Anderson, S. F. et al. 2007, astro-ph/0707.3380 

Akerman, C.J., Ellison, S.L., Pettini, M., Steidel, C.C. 2005, A&A, 
440, 499 

Anders, E., & Grevesse, N. 1989, Geochim. Cosmochim. Acta, 53, 
197 

Bohlin, R.C., Savage, B.D., & Drake, J.F. 1978, ApJ, 224, 132 
Boisse, P., Le Brun, V., Bergeron, J., & Deharveng, J.M. 1998, A&A, 
333, 841 



Calura, E, Matteucci, R, & Vladilo, G„ 2003, MNRAS, 340, 59 
Cardelli, J.A., Clayton, G.C., & Mathis, J.S. 1988, ApJ, 329, L33 
(CCM) 

Centurion, M., Bonifacio, P., Molaro, P., & Vladilo, G., 2000, ApJ, 
536, 540 

Chen, H., Kennicutt, R.C. Jr, Rauch, M. 2005. ApJ, 620, 703 
Desjacques, V, Nusser, A., & Sheth R.K. 2007, MNRAS, 374, 206 
Dessauges-Zavadsky, M., Prochaska, J.X., D'Odorico, S., Calura, F., 

Matteucci, F. 2006, A&A, 445, 93 
Efron, B. 1979, The Annals of Statistics, 7, 1 
Ellison, S.L., Hall, P.B., Lira, P. 2005, AJ, 130, 1345 
Ellison, S. L„ Yan, L„ Hook, I.M., Pettini, M., Wall, J.V, & Shaver, 

P. 2001, A&A, 379, 393 
Fall, S.M. & Pei, Y. 1989, ApJ, 337, 7 
Fall, S.M. & Pei, Y. 1993, ApJ, 402, 479 
Fontana, A., & Ballester, P. 1995, Messenger, 80, 37 
Fukugita, M., Ichikawa, T, Gunn, J.E., Doi, M., Shimasaku, K., & 

Schneider, D.P. 1996, AJ, 111, 1748 
Gordon, K.D., Clayton, G.C., Misselt, K.A., Landolt, A.U., & Wolff, 

M.J. 2003, ApJ, 594, 279 
Hampel, F.R., Ronchetti, E.M., Rousseeuw, P.J., & Stahel, W.A., 

1986, Robust Statistics. The Approach based on Influence 

Functions, (John Wiley & Sons: New York) 
Hou, J.L., Boissier, S., & Prantzos, N. 2001, A&A, 370, 23 
Jorgenson, R. A., Wolfe, A. M., Prochaska, J. X., Lu, L., Howk, J. C, 

Cooke, J., Gawiser, E., & Gelino, D. M. 2006, ApJ, 646, 730 
Junkkarinen, V. T, Cohen, R. D., Beaver, E. A., Burbidge, E. M., 

Lyons, R. W., & Madejski, G. 2004, ApJ, 614, 658 
Kulkarni, VP, Fall, S.M., Lauroesch, J.T., York, D.G., Welty, D.E., 

Khare, P., & Truran, J.W. 2005, ApJ, 618, 68 
Kulkarni, VP, York, D.G., Vladilo, G., & Welty, D.E. 2007, ApJ, 663, 

L81 

Lanzetta, K.M., Wolfe, A.M., & Tumshek, D.A. 1995, ApJ, 440, 435 
Le Brun, V, Smette, A., Surdej, J., & Claeskens, J.-F. 2000, A&A, 
363, 837 

Ledoux, C, Petitjean, P., & Srianand, R. 2003, MNRAS, 346, 209 
Linnik, Yu, V. 1961, Method of Least Squares and Principles of the 

Theory of Observations (Pergamon Press: New York) 
Lu, L., Sargent, W.L.W., Barlow, T.A., Churchill, C.W., & Vogt, S. 

1996, ApJS, 107, 475 
Menard, B., Nestor, D., Turnshek, D., et al. 2007, submitted, [astro- 

ph/0706.0898] 

Molaro, P., Bonifacio, P., Centurion, M., et al. 2000, ApJ, 541, 54 
M0ller, P., Fynbo, J.P.U., & Fall, S.M. 2004, A&A, 422, L33 
Murphy, M.T., & Liske, J. 2004, MNRAS, 354, L31 
Pei, YC. 1992, ApJ, 395, 130 

Pei, YC, Fall, S. M., & Bechtold, J. 1991, ApJ, 378, 6 
Petitjean, P., Srianand, R., & Ledoux, C. 2000, A&A, 364, L26 
Petitjean, P., Ledoux, C, Noterdaeme, P. & Srianand, R. 2006, A&A, 
456, L9 

Pettini, M., Ellison, S.L., Steidel, C.C, Shapley, A.E., & Bowen, D.V 

2000, ApJ, 532, 65 
Pettini, M., Smith, L.J., Hunstead, R.W., & King, D.L. 1994, ApJ, 426, 

79 

Pettini, M., Smith, L.J., King, D.L., & Hunstead, R.W., 1997, ApJ, 
486, 665 

Prantzos, N, & Boissier S. 2000, MNRAS, 315, 82 

Prochaska, J.X., Chen, H.-W., Wolfe, A.M., Dessauges-Zavadsky, M., 

& Bloom, J.S. 2007, ApJ, 666, 267 
Prochaska, J.X., Gawiser, E., & Wolfe, A.M. 2001, ApJ, 552, 99 
Prochaska, J.X., Hennawi, J.F., Herbert-Fort, S. 2007, ApJ, submitted 

( |astro-ph/0703594^ 
Prochaska, J.X., Herbert-Fort, S. & Wolfe, A.M. 2005, ApJ, 635, 123 



Vladilo, Prochaska, & Wolfe: color excess of QSOs with intervening DLA systems 



17 



Prochaska, J.X. & Wolfe, A.M. 1997, ApJ, 487, 73 

Prochaska, J.X. & Wolfe, A.M. 1998, ApJ, 507, 1 13 

Prochaska, J.X. & Wolfe, A.M. 1999, ApJ, 121, 369 

Prochaska, J.X. & Wolfe, A.M. 2002, ApJ, 566, 68 

Prochaska, J.X., Wolfe, A.M., Howk, J.C., Gawiser, E., Buries, S.M., 

& Cooke, J. 2007, ApJS, 171, 29 
Prochter, G.E., Prochaska, J.X., & Buries, S.M., 2006, ApJ, 639, 766 
Ostriker, J.P., & Heisler, J. 1984, ApJ, 278, 1 

Rao, S.M., Nestor, D.B., Turnshek, D.A., Lane, W.M., Monier, E.M., 

& Bergeron, J. 2003. ApJ, 595, 94 
Richards, G.T., Fan, X., Schneider, D.P., Vanden Berk, D.E., Strauss, 

M.A. et al. 2001, AJ, 121, 2308 
Schneider, D.P., Hall, P.B., Richards, G.T., Vanden Berk, D.E., 

Anderson, S.F. et al. 2005, AJ, 130, 367 
Trump, J.R., Hall, P.B., Reichard, T.A., Richards, G.T., Schneider, 

D.P., Vanden Berk, D.E., Knapp, G.R., Anderson, S.F., Fan, X., 

Brinkman, J., et al. 2006, ApJS, 165, 1 
Vanden Berk, D. E., et al. 2001, AJ, 122, 549 
Vladilo, G. 1998, ApJ, 493, 583 
Vladilo, G. 2002, A&A, 391, 407 
Vladilo, G. 2004, A&A, 421, 479 

Vladilo, G., Centurion, M., Levshakov, S.A., Peroux, C, P. Khare, 
Kulkarni, VP. & York, D.G. 2006, A&A, 454, 151 

Vladilo, G. & Peroux, C. 2005, A&A, 444, 461 

Wang, J., Hall, P.B., Ge, J., Li, A., & Schneider, D.P. 2004, ApJ, 609, 
589 

Wild, V, & Hewett, PC. 2005, MNRAS, 361, L30 
Wild, V, Hewett, P.C., & Pettini, M. 2006, MNRAS, 367, 21 1 
Wolfe, A. M., & Chen, H.-W. 2006, ApJ, 652, 981 
Wolfe, A. M., Prochaska, J. X. & Gawiser, E. 2003, ApJ, 593, 215 
Wolfe, A.M., Turnshek, D.A., Smith, H.E., Cohen, R.D., 1986, ApJS, 
61, 249 

Wolfe, A. M., Gawiser, E„ & Prochaska, J. X. 2005, ARA&A, 43, 861 
Yip, C.W., Connolly, A.J., Vanden Berk, D.E., Ma, Z., Frieman, J.A. 

et al. 2004, ApJ, 128, 2603 
York, D.G, Khare, P., Vanden Berk, D., Kulkarni, VP., Crotts, A.P.S. 

et al. 2006, MNRAS, 367, 945 



value are the same for different subsets of control quasars. The 
validity of the condition (An) = can be tested from the anal- 
ysis of quasars of the twin control samples (Section 4.1). Each 
control sample is built with unreddened QSOs and therefore its 
color deviations are An = Co - (C"). By measuring the mean 
value of these deviations we can test if (An) = 0. From Fig. 
|4] one can see that the mean values of the color excess of the 
control samples are peaked at zero, implying that the above 
condition is generally satisfied. 



Appendix A: Derivation of Eq. ©. 



By comparing the definitions E 
given in Section 3.1, we have 

A = E + A , 



= C - Co and A = C - <C> 
(A.l) 



where Ao = Co - (C"). The quantities E and Ao are physically 
independent because the color excess E depends on the phys- 
ical properties of the absorber, while Ao on the intrinsic col- 
ors of the quasars. Since E and Ao are independent variables, 
the frequency distribution of their sum equals the convolution 
product of their frequency distributions, i.e. /a = /e®/a ■ From 
this result it follows that 

(A) = (E) (A.2) 

for any arbitrary distributions and /a, provided (Ao) = 0. 
This last condition is valid for a large sample of reddened 
quasars with same properties since in this case the unknown 
unreddened color Co belongs, by construction, to the same pop- 
ulation of the unreddened colors of the control quasars C". 

If we have a large sample of reddened quasars with varying 
properties then the condition (Ao) = is still satisfied pro- 
vided the distributions of unreddened colors around their mean 



